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ABSTRACTS 


Both the total electzon content of the ionosphere (TEC) and 
the critical frequency of the P2 layer (foF2) exhibit large day- 
to-day variations during quiet and active geoaagnetic periods. It 
is of great interest to ascertain whether good correlation exists 
between TEC daily variability about the wonthly wean and foF2 
variations. With the availability of the global GPS constellation 
to provide instantaneous tiwe-delay values such a correlation say 
enable the isprovesent of HT' /•ort-ters predictions using passive 
Bonitoring of TEC. 

To detersine the correlati n, a pilot study was conducted 
using several aonths of TEC data taken in Haifa, Israel during 
1980 as %#ell as GPS tise-delay ■easureaents taken during the 
suaaer of 1992 in Jerusalea, Israel. The corresponding foF2 aea- 
sureaents were froa Cape Zevgary, Cyprus. 

The analysis showed, that for large percentages of the tine 
very good correlation exists bettreen TEC and foF2 short-tera 
variations. The correlation coefficient varies froa 0.7 - 0.8 
during winter and suaaer aonths to about 0.5 - 0.6 during equinox 
aonths. A study of the diurnal dependence of the correlation 
indicates that better correlation exists during day-tiae than 
night-tiae. There was no indication that the correlation 
coefficient is dependent on geoaagnetic activity during the 
period of this study. 

The high correlation between TEC and foF2 in the Halted data 
analysed, indicates, that real-tiae ionospheric HF prediction 
improvements are feasible when using transionospheric time-delay 
measurements. 

KEY WORDS 


Ionospheric variability; HF propagation; HF short-term pre¬ 
dictions; Total electron content (TEC); foF2; Transionospheric 
time-delay; GPS. 
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1. INTRODUCTION 

To assess the possible improvement of HF short-term predic¬ 
tions from transionospheric time-delay measurements, it is of 
great interest to ascertain whether good correlation exists 
between TEC daily variability about the monthly mean and foF2 
variability. To determine such correlation, a pilot study was 
conducted using several months of TEC data taken in Haifa, Israel 
during 1980 as well as GPS time-delay measurements taken during 
the summer of 1992 in Jerusalem, Israel. The corresponding foF2 
measurements were from Cape Zevgari, Cyprus. 

The TEC at Haifa was determined from Faraday rotation obser¬ 
vations using the signal of the Sirio satellite. The geographic 
subionospheric point corresponding to a mean ionospheric height 
of 350 km (along the signal path) is 30.3*N and 28.9*E. The 
geographic coordinates of Cape Zevgari are 34.6*N and 32.9*E. The 
measurements are thus separated by 4.3* in latitude and 4* in 
longitude. The GPS observations were taken at the National 
Physical Laboratory in Jerusalem using a Ionospheric Measurement 
System developed at the National Institute of Standards and 
Technology (NIST) at Boulder, Colorado for accurate time trans¬ 
fer. This system provides 15 minutes averages of time-delay data 
to all GPS satellites above the horizon. The equivalent TEC from 
the time-delay measurements was determined only for satellites at 
elevations larger then 30* and having a subionospheric point 
along the line of sight within 5* of the latitude of the 
ionosonde. The observed time-delays were then corrected for 
satellite biases converted into vertical TEC and averaged to 
obtain hourly values. 

2. THE BXPERHIBNTAL RESULTS 

For each month analyzed, the varability of both the TEC and 
foF2 was determined. The variability is calculated by subtracting 
the monthly average value from each hourly value and deviding by 
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the Bonthly average value. Using the hourly values of the va¬ 
riability, a cross-correlation analysis between the TEC and foF2 
was performed as well as auto-correlation for both TEC and foF2. 
Finally,the variability of the slab-thickness r, determined from 
the ratio TEC/(foF2)*, was calculated. 

The full account of the data analysis and the results obtain¬ 
ed are summerised in the paper "Assessment of foF2 Short-term 
Variations from Transionospheric Time-delay Measurements" to be 
published in the Proceedings of the International Beacon 
Satellite Symposium, which took place at Aberystwyth, Hales, from 
11-15 July 1994. The preprint of the paper is given in Appendix 
A. 


A complete set of figures for each month analyzed in this 
study are given in Appendix B. The following figures are shown 
for every month: 

a. Hourly values of the variability of foF2 and TEC. 

b. Cross-correlation function between hourly values of foF2 
and TEC. 

c. Auto-correlation functions for foF2 and TEC. 

d. Hourly values of the variability of the slab-thickness r. 

For the period June-August 1992 the results of smoothing the 
hourly values of foF2 and TEC by a 3-hour and 5-hour running 
mean, are also shown. The smoothing remarkably improved the 
correlation between foF2 and TEC, which is discussed in the 
paper in appendix A. 

3. COHCXJDSIOII 

The high cross-correlation (>0.75) coefficient for foF2 and 
TEC in the limited data presented here raises the possibility 
that real-time TEC measurements nay be used to update foF2 value 
determinations. The cross-correlation may even be higher if the 
geographic subionospheric point of the TEC measurement is closer 
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to the geographic point of the foF2 measurement, which introduces 
an error, in addition to the possible inherent measurements un¬ 
certainties. TEC measurements utilizing satellite emitted 
signals are passive in nature and do not burden the 
electromagnetic spectrum. In addition, the availability of the 
global GPS constellation to provide instantaneous time-delay, or 
equivalently TEC, values could provide an instantaneous updating 
of foF2 models on a global basis as well as on a regional basis. 
Such capability is important for HF communication along short, 
medium and long range paths. 
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ASSSSSIOBST OF fgFj SK«T>TEIUf VAKIATZONS FROM 
TRANSIONOSFBERIC TZMB-OZXAY MERSUlUEMEirrS 

Z. Kouminer^ and H. Soicher^ 

1. Asher Space Research Institute, Technion - Israel Institute of Technclog-y, 

Haifa 32000, Israel 

Tax: *972 4 230956, e-mail: as2wih9vjnsa.technion.ac.il 
2. Space and Terrestrial Communications Directorate, 

US Army Communications-Electronics Command, Fort Monmouth, 

NJ 07703-5203, USA 

Fax 1(908)532-0456, e-mail: Soicher%DOIM69Monmouth - emh3.arrtiy.mil 

1. INTRODUCTION 

To assess the possible improvement of HF short-term predictions from 
transionospheric time-delay measurements, it is of great interest to ascertain 
wnether good correlation exists between TEC daily variability about the 
monthly mean and foF 2 variaOoility. To determine such correlation, a pilot 
study was conducted using several months of TEC data taken in Haifa, Israel 
during 1980 (Soicher et al., 1982) as well as GPS time-delay measurements 
ta)cen during the summer of 1992 in Jerusalem, Israel. The corresponding foF 2 
measurements were from Cape Zevgari, Cyprus. 

The TEC at Haifa was determined from Faraday rotation observations using the 
signal of the Sirio satellite. The geographic subionospheric point 
corresponding to a mean ionospheric height of 350km (along the signal path) is 
30.3°N and 28.9°E. The geographic coordinates of Cape Zevgari are 34.6°N and 
32.9°E. The measurements are thus separated by 4.3® in latitude and 4° in 
longitude. The GPS observations were taken at the National Physical 
Laboratory in Jerusalem using a Ionospheric Measurement System developed at 
the National Institute of Standards and Technology (NIST) at Boulder, Colorado 
for accurate time transfer. This system provides 15 minutes averages of time- 
delay data to all GPS satellites above the horizon. The equivalent TEC from 
the time-delay measurements was determined only for satellites at elevations 
larger than 30° and having a subionospheric point along the line of sight 
within 5° of the latitude and longitude of the ionosonde. The observed time 
delays were then corrected for satellite biases using the JPL table of 
corrections (Wilson and Mannucci, 1993), converted into vertical TEC 
(Klobuchar, 1987) and averaged to obtain hourly values. 

2. EXPERIMENTAL RESULTS 

Figure 1 shows the hourly values of the variability of foF 2 and TEC for the 
period 1-16 August 1992. That period is near the maximum of solar cycle 22 
(monthly mean sunspot number 103). The variability is determined by 
subtracting the monthly average value from each hourly value and dividing by 
the monthly average value. The results of cross correlation analysis on the 
and TEC depicted in Figure 1 are shown in Figure 2. A maximum cross- 
correlation coefficient of 0.7 occurs at zero time lag while the correlation 
reduces very quickly with time lag. It is thus shown chat the correlation 
between foF 2 and TEC is very good. Auco-correlacion results for both fof"2 
TEC for the same time period (Figure 3) shows that the auto-correlation is 
similar in character to the cross correlation function and drops to 0.5 after 
a time lag of 3-4 hours. This is consistent with results obtained elsewhere 
indicating that short-term predictions based on real-time observations can 
only be useful for a few hours. 
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It can be seen from Figure 1 that both the variations in foF 2 and T£C are 
rather noisy on an hour-to-hour basis which may be caused by measurements and 
data reduction errors rather than by physical phenomena. The results or 
smoothing the variations in TEC and foP'2 * 3-hour running mean is shown in 
Figure 4. It can be seen that the correlation is remarlcably improved 
(correlation coefficient of 0.8). This shows that indeed the data are noisy 
and that the correlation between TEC and foF 2 is actually better chan what the 
raw data indicate. 

The seasonal dependence of the correlation coefficient between the variability 
in TEC and fQF 2 is shown in Table 1. It can be seen that the correlation 
coefficient is better chan 0.7 in winter and summer. However, during equinox, 
especially in April and May 1980, the correlation drops to about 0.5. The 
reason for the low correlation is the large post sunset enhancements in TEC 
which were not seen in ^ 0 ^ 2 - These post sunset enhancements, observed mainly 
in equinox months, are attributed to electron fluxes arriving from the 
equatorial regions along the magnetic lines of force. They are latitude- 
dependent (Soicher ec al., 1984) and thus did not affect the foF2 
observations. 

The diurnal dependence of the correlation is shown in Table 2. It can be seen 
that higher correlation occurs during day time cheuo night time. Between 10-18 
hours local time, the correlation coefficient is about 0.7; while between 22- 
06 hours local time, the cor ^ation drops to between 0.55-0.6. 


TABLS-l TAflL£-2 


SEASON 

MONTH 

CORRELATION 

COEFFICIENT 

HOURS LT 

CORRELATION 

COEFFICIENT 

Winter 

Jan 

80 

0.74 

02-06 

0.55 


Feb 

80 

0.73 

06-10 

0.63 

Equinox 

Mar 

80 

0.66 

10-14 

0.70 


Apr 

80 

0.50 

14-18 

0.71 


Mav 

80 

0.47 

18-22 

0.60 

Summer 

Jul 

81 

0.73 

22-02 

0.61 


Jul 

92 

0.77 




Aug 

92 

0.70 




The ionosphere is )cnown to vary substantially with magnetic activity (Goodman, 
1992) . To ascertain whether magnetic activity has any impact on the cross 
correlation of fof ‘2 TEC, two time periods - one quiet (18-26 January 1980) 
and one active (14-19 Febn;ary 1980) - were examined. The former with index 
Ap<6 is depicted in Figure 5 with maximum correlation coefficient of 0.823 and 
the latter with index Ap<40 is depicted in Figure 6 with maximum correlation 
coefficient 0.774. Thus it appears that magnetic changes do not have a mar)ced 
impact on the correlation of the two ionospheric parameters. 

The slab thic)<ness T is determined from the ratio TEC/(fQF 2 )^- It is expected 
that the variability of T will be smaller than the variability of either the 
TEC or fof2 because of the good correlation between these two ionospheric 
parameters. Thus, global models of the slab thic)aiess updated with real-time 
measurements of TEC obtained with the GPS networ)t, might give improved values 
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of foF'2- However, the calculated variability in x is of the same order as the 
corresponding variability of TEC or fo^2- example, for the period 11-29 

February 1980, the standard deviation of the hourly variabi''ity of x is 16% in 
comparison to that of TEC which is 22%. 

3. CONCLUSIONS 

The high cross-correlation (>0.75) coefficient for foF 2 and TEC in the limited 
data presented in this paper raises the possibility that real time TEC 
measurements may be used to update foF ‘2 value determinations. The cross¬ 
correlation may even be higher if the geographic subionospheric point of the 
TEC measurement is closer to the geographic point of the fof ‘2 measurement, 
which introduces -n error, in addition to the possible inherent measurement 
uncertainties. 1 measurements utiliring satellite emitted signals are 

passive in nature and do not burden the electromagnetic spectrum. In 
addition, the availability of the global GPS constellation to provide 
instantaneous time-delay, or equivalently TEC, values could provide an 
instantaneous updating of foF 2 models on a glroal basis as well as on a 
regional basis. Such capability is important for HF communication along 
short, medium and long range paths. 
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Soicher, H., Z. Houminer and A. Shuval, Total Electron content in the Middle 
East, Radio Science, ^7, 1623-1631, 1982. 

Soicher, H., J. A. Klobuchar and P. H. Doherty, Spatial Variability of Total 
Electron Content in the Eastern Mediterranean Region, Radio Science, 

757-754, 1984. 

Wilson, 3. D. and A. J. Mar.nucci, Instrumental Biases in Ionospheric 
.Measurements Derived from GPS Data, Proceedings of the Institute of Navigation 
GPS-93, Salt Lake City, Utah, September 1993. 

5. FIGURE CAPTIONS 

1. Hourly values of the variability of f^F 2 (at Cyprus) and TEC (at Israel) 
for the time period 1-16 August 1992. 

2. Crosscorrelation function for the foF2 values of Figure 1. 

3. Autocorrelation functions for the foF 2 and TEC values of Figure 1. 

4. Smoothed hourly values of the variability values of Figure 1. 

5. The variability (in percent) of fQF 2 and TEC for the geomagnetically quiet 
period of 18-26 January 1980. 

6. The variability (in percent) of foF 2 and TEC for the geomagnetically 
active period of 14-19 February 1980. 
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APPENDIX B 


Experinental results for the period January - May 1980, 
July 1981 and June - August 1992. 

The following are the figure captions: 


Fig B-1 
Fig B-2 
Fig B-3 
Fig B-4 
Fig B-5 
Fig B-6 
Fig B-7 
Fig B-8 
Fig B-9 
Fig B-10 
Fig B-11 
Fig B-12 
Fig B-13 
Fig B-14 
Fig B-15 


Hourly values of the variability of foF2 and TEC 
for the period 8>31 January 1980. 
Cross-correlation function between foF2 and TEC 
for the period 8-31 January 1980. 
Auto-correlation functions for foF2 and TEC for 
the period 8-31 January 1980. 

Hourly values of the varieibility of r for the 
period 8-31 January 1980. 

Hourly values of the variability of foF2 and TEC 
for the period 11-29 February 1980. 
Cross-correlation function between foF2 and TEC 
for the period lJL-29 February 1980. 
Auto-correlation functions for foF2 and TEC for 
the period 11-29 February 1980. 

Hourly values of the variability of t for the 
period 11-29 February 1980. 

Hourly values of the variability of foF2 and TEC 
for the period 1-31 March 1980. 

Cross-correlation function between foF2 and TEC 
for the period 1-31 March 1980. 

Auto-correlation functions for foF2 and TEC for 
the period 1-31 March 1980. 

Hourly values of the variability of t for the = 
period 1-31 March 1980. 

Hourly values of the variability of foF2 and TEC 
for the period 1-30 April 1980. 

Cross-correlation function between foF2 and TEC 
for the period 1-30 April 1980. 

Auto-correlation functions for foF2 and TEC for 
the period 1-30 April 1980. 
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Pi9 B-16 
Fig B-17 
Pig B-18 
Pig B-19 
Fig B-20 
Fig B-21 
Fig B-22 
Fig B-23 
Fig B-24 
Fig B-25 
Fig B-26 
Fig B-27 
Fig B-28 
Fig B-29 
Fig B-30 
Fig B-31 
Fig B-32 
Fig B-33 


Hourly vuluos of tho variability of r for tha 
pariod 1>30 April 1980. 

Hcmrly valuas of tha variability of foF2 and TEC 
for tha pariod 1-31 May 1980. 

Crosa-corralation function batwaan foF2 and TEC 
for tha pariod 1-31 Nay 1980. 

Auto-corralation functiona for foF2 and TEC for 
tha pariod 1-31 May 1980. 

Hourly valuaa of tha variability of r for tha 
pariod 1-31 Nay 1980. 

Hourly valuas of tha variability of foF2 and TEC 
for tha pariod 2-31 July 1981. 

Cross-correlation function between foF2 and TEC 
for the period 2-31 July 1981. 

Auto-correlation functions for foF2 and TEC for 
the period 2-31 Ju^y 1981. 

Hourly values of the varied>ility of r for the 
period 2-31 July 1981. 

Hourly values of the variability of foF2 and TEC 
for the period 1-30 June 1992. 

Smoothed hourly value by a 3 hour running mean 
for the period 1-30 June 1992. 

Smoothed hourly value by a 5 hour running mean 
for the period 1-30 June 1992. 

Cross-correlation function between foF2 and TEC 
for the period 1-30 June 1992. 

Auto-correlation functions for foF2 and TEC for 
the period 1-30 June 1992. 

Hourly values of the variability of r for the 
period 1-30 June 1992. 

Hourly values of the variability of foF2 and TEC 
for the period 2-31 July 1992. 

Smoothed hourly value by a 3 hour running mean 
for the period 2-31 July 1992. 

Smoothed hourly value by a 5 hour running mean 
for the period 2-31 July 1992. 
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Pig B>34 
Pig B-35 
Pig B-36 
Pig B-37 
Pig B-38 
Fig B-39 
Pig B-40 
Fig B-41 
Fig B-42 


Cross-corrslation function botwoon foP2 and TEC 
for tho poriod 2-31 July 1992. 

Auto-corrolation functicms for foP2 and TEC for 
tha pariod 2-31 July 1992. 

Hourly values of tha variability of r for tha 
pariod 2-31 July 1992. 

Hourly values of tha variability of foP2 and TEC 
for tha pariod 1-16 August 1992. 

Saoothad hourly value by a 3 hour running saan 
for tha pariod 1-16 August 1992. 

Saoothad hourly value by a 5 hour running aean 
for tha pariod 1-16 August 1992. 
Cross-correlation function between foF2 and TEC 
for tha pariod 1-16 August 1992. 

Auto-corralatlon functions for foF2 and TEC for 
tha pariod 1-16 August 1992. 

Hoiurly values of tha variability of r for the 
pariod 1-16 August 1992. 
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8-31 January 1980 



TIME DELAY [hours] 


Fig B-5 
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11-29 February 1980 

C = 0.733 
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11-29 February 1980 



Fig B-6 
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11-29 February 1980 



TIME DELAY [hours] 


Fig B-7 
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11-29 February 1980 
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Fig B-8 
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Fig B-11 
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1-30 April 1980 

C=0.505 



TIME [hours] 


Fig B-l.> 




















Fig B-14 
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1-30 April 1980 



TIME DELAY [hours] 


Fig B-15 
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1-31 May 1980 

C = 0.471 



TIME [hours] 


Fig B-17 
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1-31 May 1980 



Fig B-IS 
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1-31 May 1980 



TIME DELAY [hours] 


Fig B-19 
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1-31 May 1980 



TIME [hours] 


Fig B-20 
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2-31 July 1981 



TIME DELAY [hours] 


Fig B-2 
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1 - 30 June 1992 
3 point smoothing C = 0.762 



TIME [hours] 


Fig B-26 
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2-31 July 1992 

C-0.778 



Fig 
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2-31 July 1992 
3 point smoothing C = 0.835 
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2-31 July 1992 

5 point smoothing C = 0.859 



Fig 
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2 - 31 July 1992 



TIME DELAY [hours] 


Fig B-55 
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1-16 August 1992 



Fig B-41 






































